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Abstract. We present model calculations to the observed 
spectral variability in the 7-ray spectrum of the ultralu- 
minous blazar-type quasar PKS 0528+134. We argue that 
the observed appearance of a spectral break between the 
COMPTEL and the EGRET energy range can plausibly 
be explained by a variation of the Doppler beaming fac- 
tor in the framework of a relativistic jet model for AGNs, 
combining the ERG and the SSG model. This may be an 
instructive example for the different dependence of these 
two spectral components on the Doppler factor, as pre- 
dicted in earlier work. 
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1. Introduction 

Up to now, more than 50 blazar-type AGN have been de- 
tected by EGRET as emitcrs of high-energy 7-rays above 
100 MeV (Mattox et al. 1997). These sources are identi- 
fied with flat-spectrum (a > -0.5) radio sources classified 
as QSOs and BL Lac objects. A large fraction of these 
blazars exhibits variability at 7-ray energies on time scales 
of days to months (Mukherjee et al. 1997). At radio wave- 
lengths, all blazars can be recognized as bright, compact 
sources with a flat synchrotron spectrum emanating from 
outflowing plasma jets that are nearly aligned with our 
line of sight. Relativistic beaming is required in the ob- 
jects in view of the luminosity and variability time scales 
(Dermer and Gehrels 1995), in accord with VLBI obser- 
vations indicating that superluminal motion is a common 
feature in this class of AGN (e. g., Pohl et al. 1995, Barthel 
et al. 1995, Krichbaum et al. 1995). 

One of the most luminous examples of this class of 
obiects is the auasar PKS 0528-1-134. It is known as a 
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VLBI images show a one-sided core-jet structure of ^ 5 
mas length and superluminal motion (Pohl et al. 1995). 
An optical spectrum, obtained in November 1991 with the 
ESO 3.6 m telescope at La Silla / Ghile, revealed a red- 
shift of z = 2.07 (Hunter et al. 1993). At X-rays, the source 
was first detected with the Einstein observatory (Bregman 
et al. 1985). At 7-ray energies, PKS 0528+134 is detected 
by EGRET, COMPTEL and OSSE aboard the Compton 
Gamma-ray Observatory (GGRO) (Mukherjee et al. 1996, 
Collmar et al. 1997, McNaron-Brown et al. 1995), showing 
a time-variable 7-ray flux with intensity variations up to 
a factor of ^20 (Mukherjee et al. 1996). The broadband 
spectrum of PKS 0528+134 reveals that during 7-ray high 
states its total energy output is dominated by its 7-ray 
emission. Under the assumption of isotropic emission, the 
total luminosity of the object during 7-ray high states 
ergs s~^ (throughout this paper, we 
75 km s"^ Mpc"^ and qo 0.5). 

The strongest EGRET blazar detections can be char- 
acterized by a single power-law spectrum with differential 
photon spectral indices 1.5 < s < 2.7 (Thompson et al. 
1995). In some cases, such as PKS 0528+134 (Mukherjee 
et al. 1996), the EGRET spectrum (30 MeV - -10 GeV) 
seems to harden during 7-ray flares. A recent analysis of 
the first 3.5 years of COMPTEL observations of PKS 0528 
+134 (Collmar et al. 1997) indicated a correlation between 
its 7-ray flux and the occurence of a spectral break at 
MeV energies. During states of high 7-ray intensity the 
combined EGRET and COMPTEL spectra (-1 MeV to 
~10 GeV) are characterized by a spectral break between 
the EGRET and the COMPTEL energy range, while dur- 
ing 7-ray low states the 7-ray spectrum from MeV to GeV 
energies is consistent with a single power-law. 

Measurements of broadband spectral variability in 7- 
ray blazars can put the most severe constraints on mod- 
els for these objects. By now, it is widely accepted that 
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positrons) in a relativistic jet directed at a small angle 
with respect to our line of sight (e. g. Schlickeiser 1996). 

A major subject of debate is the source of the soft pho- 
tons which are upscattered to 7-ray energies. These can 
be synchrotron photons produced by the jet itself (the 
SSC model, e. g. Maraschi et al., 1992, Bloom & Marscher 
1996), accretion disk radiation entering the jet directly 
(the ERC model, Dermer, Schlickeiser & Mastichiadis 
1992, Dermer & Schlickeiser 1993) or after being rescat- 
tered by the broad line regions (Sikora, Begelman & Rees 
1994), or jet radiation being reflected by a torus of cold 
material and reentering the jet (GhiscUini & Madau 1996). 

Detailed investigations of the ERC model and of the 
SSC model have recently been presented by Dermer, Stur- 
ner & Schlickeiser (1997) and by Bloom & Marscher 
(1996), respectively. These reveal that sharp spectral 
breaks in the 7-ray spectra of blazars at several MeV fa- 
vor the ERC mechanism as a possible explanation because 
such breaks cannot be produced by a homogeneous SSC 
model. In contrast, a smooth transition from a hard spec- 
trum at X-ray and soft 7-ray energies to a softer EGRET 
spectrum seems indicative of SSC being the dominant ra- 
diation mechanism. Nevertheless, a pure version of one 
of these models seems to be an unrealistic simplification 
since both of these radiation mechanisms are plausibly at 
work in relativistic AGN jets. A synthesis of both mod- 
els has been investigated in detail by Bottcher, Mause & 
Schlickeiser (1997). 

Dermer et al. (1997) have shown that the ERC radi- 
ation exhibits a stronger dependence {Ferc{^) k D^^^) 
on the Doppler factor D = (F [1 — fir cosOobs])^^ than 
the SSC component {Fssc{^) oc d(5+«)/2). Here, T = 
(1 — /3r)^^''^ is the bulk Lorentz factor of the ultrarela- 
tivistic jet plasma and s is the spectral index of the initial 
particle distribution function. The kind of spectral vari- 
ability as observed in PKS 0-528-1-134 can therefore easily 
be explained by a change in the Doppler factor (i. e. a 
change in the bulk Lorentz factor or in the angle Oobs be- 
tween the direction of motion and the line of sight) in the 
framework of a combined SSC/ERC model. According to 
this idea, we suggest that in the case of PKS 0528+134 
the SSC radiation dominates the 7-ray spectrum in the 
7-ray low state (with no spectral break in the COMPTEL 
regime required), while in the 7-ray high state an addi- 
tional ERC component becomes dominant in the 7-ray 
spectrum at energies around 10 MeV. In this Letter, we 
present model calculations demonstrating the viability of 
this idea. 

2. The model 

To produce spectral fits to PKS 0528+134 we use the com- 
bined ERC/SSC model described in detail by Bottcher et 



black hole of M = 5-lO^°M0. The luminosity of the accre- 
tion disk is assumed to be io = 5 • lO^^ergs"^. The blob 

moves with a constant bulk Lorentz factor P. The parti- 
cles inside the blob are isotropically distributed according 
to a power-law 71(7) oc 7"'* for 71 < 7 < 72 in the rest 
frame of the blob. We assume reacceleration in the initial 
phase to be inefficient in comparison to radiative cooling 
of the pairs and follow the self-consistent evolution of the 
pair distribution as the blob moves out. taking into ac- 
count synchrotron radiation, SSC scattering to arbitrarily 
high order, ERC scattering, 7-7 absorption and pair pro- 
duction. Comparison to observations is made on the basis 
of the time-averaged flux of a blob since the inferred cool- 
ing timescalc of the material in the blob (in the observer's 
frame) is shorter than the time resolution of COMPTEL 
and EGRET. 

In our model we assume the ejection of individual plas- 
ma blobs with different bulk Lorentz factors and slight 
differences in the internal energy distributions of the ma- 
terial in the blob. These are described by the following 
parameters: 71^2 (the low- and high-energy cutoffs in the 
Lorentz factors of the electrons), s (the power-law index 
of the electon distribution), B (the magnetic fleld), F (the 
bulk Lorentz factor of the blob), Oobs (the angle between 
the line of sight and the jet axis), z (the starting height 
of the center of the blob above the accretion disk). 




v[Hz] 



Fig. 1. Fit to the (non simultaneous) broadband spectrum of 

PKS 0528+134. Parameters: 71 = 2 • 10^ 72 = 9 ■ 10^, s = 3.0, 
n = 60cm~^ Rb = A ■ 10^** cm, B = 1.3 G, z = 4 • 10"^ pc, 
Lo = 5 • 10'"'crgs-\ r = 20, = 2.5°. Dashed: ERC, 
dot-dashed: SSC, long dashed: synchrotron, dotted: accretion 
disk, solid: total 



3. Spectral fits 
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the broadband spectrum as well as the peculiar spec- 
tral variability of this object. As a first step, we produce 
a fit to the (non simultaneous) broadband spectrum of 
PKS 0528+134 which is shown in Fig. 1. The radio data 
are taken from Reich et al. (1993). The optical data point 
is still imder debate because PKS 0528+134 is located be- 
hind a strongly absorbing cloud of interstellar material, 
and the correct dereddening factor is unclear (Wagner 
1997, private communication). Here, we use the same opti- 
cal point as used by Mukherjee et al. (1996). The ROSAT 
point is taken from Zhang et al. (1994), and the OSSE data 
point is from McNaron-Brown et al. (1995). The contem- 
poraneous COMPTEL and EGRET data from Viewing 
Period (VP) (April 1991) are taken from Collmar et al. 
(1997). 

Even though the broadband spectrum shown in Fig. 
1 is not measured simultaneously, it gives us an idea of 
the parameters to be used in the detailed fits to the 7- 
ray spectra in the high and low states of PKS 0528+134, 
as presented by Collmar et al. (1997). In the model de- 
scribed above, it is assumed that the radio spectrum (be- 
low ~ 100 GHz) is predominantly produced in the more 
distant parts of the jet, because the 7-ray emitting plasma 
is optically thick to synchrotron self-absorption at radio 
frequencies. The OSSE data point is not reproduced by 
our fit because it has not been measured contemporane- 
ously with the 7-ray spectrum and docs therefore not cor- 
respond to the 7-ray high state. (Figure 2, however, illus- 
trates that our spectral fit to the 7-ray low state is in good 
agreement with the OSSE spectrum at that time.) 

Fig. 1 demonstrates that the spectral break between 
the COMPTEL and the EGRET energy range can well 
be explained by an ERG component dominating the spec- 
trum at energies around 10 MeV in the 7-ray high state. 
The strong spectral break being located at ~ 10 MeV con- 
strains the input soft photons to the component causing 
this break to a relatively low energy. Given the required 
total luminosity of the disk photon field, this leads us to 
the conclusion of an accretion disk radiates at only w 1 % 
of the Eddington luminosity for a central black-hole mass 
of ~ 5 • 10^° Mq, which is consistent with dynamical stud- 
ies of the central region of PKS 0528+134. The size of the 
emitting region of Rb ~ 4 • 10^^ cm is consistent with the 
timescale of ~ 2 days for intensity variations of order 100 
% as found recently by Wagner et al. (1997). 

The low state, as shown in Fig. 2, could well be repro- 
duced by a pure SSC model. For our fit to this state we 
used basically the same set of parameters as for the high 
state (Fig. 1), but a lower bulk Lorentz factor (F = 7 in 
the low state versus F = 20 in the high state). The differ- 
ent dependence of the SSC and the ERG component on 
the Doppler factor (Dermer et al. 1997) then causes the 
observed spectral variability in the 7-ray regime. Since we 
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Fig. 2. Fit to the 7-ray (OSSE: triangles, COMPTEL: circles, 
EGRET: squares) spectrum of PKS 0528+134 in its low state 
(combined VP 321/322). Parameters are the same as for Fig. 
1, except r = 7, s = 2.8. Dashed: ERG, dot-dashed: SSG, solid: 
total 

equally good as the one shown in Fig. 1) does not yield any 
physically relevant information. In all figures the SSG and 
ERG components are plotted without taking 7-7 absorp- 
tion into account. In the calculation of the total spectrum 
as well as for the simulation of the evolution of the jet 
material, 7-7 absorption and pair production are included 
self-consistently. This is the reason for the (unabsorbed) 
SSC contribution being higher than the total escaping flux 
at high energies in the figures. 
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Fig. 3. Fit to the 7-ray (GOMPTEL and EGRET) spectrum 
of PKS 0528+134 in its high state in VP 1. Parameters are the 
same as for Fig. 1, except 71 = 2.4 • lO"*, 72 = 8 ■ 10'', s = 3.1, 
r = 22. Dashed: ERG, dot-dashed: SSG, solid: total 
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factor as compared to the fit shown in Fig. 1, lead to a 
good fit to the 7-ray high states observed in viewing pe- 
riod 1 , while the state showing the highest flux in the his- 
tory of these observations, found in VP 213, seems indeed 
to be exceptional and requires a very high bulk Lorentz 
factor and a significantly harder particle spectrum than 
the other states. Its 7-ray spectrum could equally well be 
reproduced with a pure ERC model, thus confirming the 
trend that the more luminous the object is in 7-rays the 
more dominant is the ERC component. The relevant pa- 
rameters of the different fits are given in Table 1. 
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Fig. 4. Fit to the 7-ray (COMPTEL and EGRET) spectrum 
of PKS 0528+134 in its high state in VP 213. Parameters are 
the same as for Fig. 1, except 72 = 10^, s = 2.4, V = 25. 
Dashed: ERC, dot-dashed: SSC, soUd: total 



Table 1. Interesting model parameters for the different obser- 
vational periods. -Fj_jo4 is the integrated flux from 1 MeV 
to 10 GeV of our model fits in units of keV cm~^ s~^. The dif- 
ferent 7-ray fluxes reflect mainly the different Lorentz factors 
r of the individual blobs. 



Parameter 


VP 321/322 


VP 


VP 1 


VP 213 


71 


2 • 10^ 


2 • 10^ 


2.4 ■ 10^ 


2 • 10^ 


72 


9-10* 


9-10'' 


8- lO"* 


10^ 


s 


2.8 


3.0 


3.1 


2.4 


r 


7 


20 


22 


25 


-f'l-lO*' 


0.6 


1.1 


1.1 


3.3 



4. Summary and conclusions 



exhibits a peculiar spectral behavior between its 7-ray low 
and high states, indicative of an additional spectral com- 
ponent showing up around several McV in the high state. 

Our model fits demonstrate that a combined ERC / 
SSC model is very well suited to reproduce this behavior, 
based on the different dependence of the ERC and the SSC 
component on the Doppler boosting factor. In the 7-ray 
low state, the 7-ray spectrum can basically be reproduced 
as a pure SSC spectrum, while in the 7-ray high state, 
due to an increase in the bulk Lorentz factor of the jet, 
the ERC component becomes dominant over a significant 
part of the 7-ray spectrum, producing a spectral break 
around several MeV. 

The spectral variability in PKS 0528+134 appears to 
be an interesting manifestation of the different Doppler 
boosting patterns of the ERC and the SSC radiation, as 
predicted earlier. A similar variation of the Doppler factor 
can, of course, also result from a change in the viewing 
angle. For example, the reduction of the beaming factor as 
appropriate to explain the spectrum of VP 321/322 (Fig. 
2) can equally result from assuming a viewing angle of 
4.7° (instead of 2.5° as used for our simulations) when 
leaving the bulk Lorentz factor constant at F = 20. The 
diagnostic on the Doppler factor presented in this Letter 
cannot distinguish between these two effects. 
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